Radiation back reaction on moving branes 
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This paper addresses the radiation back reaction problem for cosmological branes. A general 
framework is provided in which results are given for the radiation reaction with massles and mas- 
sive scalar fields with flat extra dimensions and massless conformal fields in anti-de Sitter extra 
dimensions. For massless scalar field radiation the back reaction terms in the equation of motion 
are non-analytic. The interpretation of the radiation reaction terms is discussed and the equations 
of motion solved in simple cases. Nucleosynthesis bounds on dark radiation give a lower bound on 
the string vacuum energy scale of V 'At m p , where At is the tensor perturbation amplitude in the 
cosmic microwave background. 



I. INTRODUCTION 



Attempts to construct a unified theory of gravity and the other fundamental forces using superstrings have provided 
a new impetus for the study of higher dimensions and branes. In certain low energy limits, superstring theory reduces 
, to a brane universe, or universes, where ordinary matter is confined to surfaces embedded in higher dimensions pj. 
The cosmological evolution of these models can be different from the standard cosmological scenario L3- Li LJl and their 
study is therefore worthwhile. 

The matter fields on the brane interact with higher dimensional fields, which include the graviton and any particles 
associated with the graviton by supersymmetry. Cosmological equations, which describe the evolution from the 
point of view of the brane, can be obtained by decomposing the higher dimensional Einstein equations and applying 
boundary or junction conditions for the brane 0. These equations reduce to the standard four dimensional equations 
tyQ . for energies small compared to the vacuum energy scale of the brane. 

Our interest here is mainly in the effects of radiation generated by the movement of the brane. We shall focus 
our attention on the radiation back reaction force on the brane. The radiation is closely related to the radiation due 
to a moving mirror p|. The mirror, like the brane, disturbs the vacuum fluctuations of the radiation field through 
the effect of the boundary condition at the surface of the mirror. The vacuum fluctuations react back on the mirror 
■ though the pressure component in the stress energy tensor. 

It may be useful to recall some of the old results on moving mirrors with massless scalar field radiation. The 
pressure force is proportional to the 3'rd time derivative of the mirror's position in two dimensions [j| and the 5'th 
1 time derivative in four dimensions Q. There is no radiation reaction force on a uniformly accelerating mirror. A 
^ , rather perplexing result is that the pressure component of the stress energy tensor diverges as the distance from the 
mirror shrinks to zero . It is believed that the correct pressure on the mirror is obtained by dropping the divergent 
terms. Some support for this view is provided by the fact that the total work done by the divergent terms is zero, at 
jf^ ■ least when the mirror motion has constant velocity asymptotically. 

The moving mirror might also be placed in a heat bath with temperature T. In this case, for a mirror moving with 
a speed v, the thermal contribution to the pressure force is proportional to T 2 v in two dimensions and T 4 v in four 
dimensions 8] . These agree with the forces which would be expected according to elementary kinetic theory applied 
to a mirror and a gas of photons. 

The equation of motion for the mirror can be modified to take these radiation reaction forces into account. The 
modified equation of motion is analogous to the Abraham-Lorentz equation of motion for an electron with radiation 
damping [9| , which offers a guide to the interpretation of the radiation reaction effects for the moving mirror and the 
moving brane. The Abraham-Lorentz equation for the electron's momentum is 

where F is an external force and r = 2e 2 /3mc 3 . A feature of this equation is the existence of 'runaway' solutions 
p oc exp(t/r), which are rejected on the grounds that the radiation reaction term should be only a small correction to 
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the equation of motion. The runaway solutions can be eliminated by replacing the equation by an integrodifferential 
equation, but this introduces another feature, 'preacceleration', where the electron begins to respond a short time 
before the external force is applied. We shall see that these features can be found in the brane system. 

Recent work related to radiating branes includes ^3 > which examined the radiation due to a uniformly accelerated 
brane. The authors did not consider the radiation reaction, but it would be expected to vanish for this type of motion. 
The brane can also radiate though interactions of the type X + X — ► Y, where the X's are particles of matter and 
Y is a bulk graviton. The effect of these gravitons on the brane persists, and resembles ordinary radiation. In some 
very specific models of brane cosmology, the gravitons can make a small contribution to the expansion rate during 
the nucleosynthesis era. This has been investigated by a number of authors 0, 0, ^| . A similar type of reaction, 
where X is the inflaton and Y a bulk scalar, could also change the nature of reheating after a period of inflation 0] . 

The next section contains some general results for the radiation from moving branes and mirrors. We take some 
care with the perturbation theory because we would like to ensure that the calculations give the vacuum expectation 
value (in|T a fc|in) of the stress energy tensor rather than (out|T a f,|in). It is also important to separate the radiation 
reaction calculation from effective action and casimir energy calculations. From the point of view of the brane and its 
cosmological evolution, the radiation back reaction effect is a dissipative phenomenon and cannot be derived from a 
reduced effective action which only depends on the degrees of freedom attached to the brane. However, we shall see 
that the Schwinger-Keldeysh formalism gives a reduced effective action which can be used to analyse both dissipative 
and non-dissipative effects. 



II. RADIATION BACK REACTION 



A. Backreaction equations 



Our approach to the back reaction on a brane moving in one extra dimension is closely based on the back reaction 
problem for a moving mirror In both cases the radiation is produced because the radiation satisfies boundary 
conditions on a moving boundary. In moving mirror problems the back reaction force is determined by the pressure 
components of the expectation value of the stress energy tensor. 

The back reaction of this radiation on the motion of a brane in five dimensions can be obtained by reducing the five 
dimensional Einstein equations in the manner described by Shiromizu, Maeda and Sasaki Q. This reduction assumes 
a reflection symmetry between to two sides of the brane. For the 5-dimensional Einstein equations we have 

Gab = K^Tab (2) 

For simplicity, we shall consider just a 5-dimensional cosmological constant A 5 , a brane vacuum energy A and 5- 
dimensional radiation, 

Tab = 2 A 5 g ah - \h ab 8(E) + (m|T£ 6 |m) (3) 

where h a b is the metric induced on the brane by the 5-dimensional metric g a b- Using the Gauss-Codacci equations, 
Shiromizu, Maeda and Sasaki show that the Einstein tensor ^G^ satisfies 

(4) G F + AV = + {T r nn - -T r ) h^\m) E„ v (4) 

where E^ v is a projection of the Weyl tensor, C n ^ nv and therefore trace-free. The 4-dimensional cosmological constant 
term and Newton's constant are related to their 5-dimensional counterparts by 

A = ^(6A 5 + «SA") (5) 

G = t < 6 > 

In order to recover the correct low energy limit for the single brane we require A > and, for a small or vanishing 
cosmological constant, A5 < 0. 

The E^ v term can be eliminated by taking the trace of l@J), giving an equation for the Ricci scalar of the brane, 

( 4 )i? = 4A-2^( in |T;jm) (7) 



In a cosmological context, this is the equation of motion of the scale factor. The expectation value on the right 
includes the radiation damping effect on the expansion rate of the universe. In general, the expectation value depends 
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on both the five dimensional metric and the brane motion. However, in this paper we shall assume that the back 
reaction of the radiation on the five dimensional metric can be neglected in the quantum calculation. 

For a spatially flat, homogeneous universe with scale factor a(t) and expansion rate Hit), the cosmological evolution 
equation 10 becomes 

6H + 12H 2 = 4A - 2«|(in|7£ n |in). (8) 

The initial condition as t — > — oo is provided by the Friedman equation 3H 2 = A. The Friedman equation is modified 
at all other times by the Weyl term, 

E 00 = ^ (in] (T^ - r;j|in) + A - 3H 2 , (9) 

which can be interpreted as energy lost by the brane. Note, however, that the quantum term can be negative. 

Typically, the radiation reaction (in|T^ n |in) will only be significant for very early times, but it leaves behind a 
residual 'dark radiation' term 0] obtained by integrating JSJl, 

2k- 2 r°° 

E 00 ~-t Hit'Mt'fWTU^Wdt' (10) 
a J-oo 

In principle, this term behaves like a contribution to the radiation energy density, but we have not, so far, considered 
the effects of ordinary matter. If, for example, matter is generated by reheating after a period of inflation, then the 
dark radiation produced during or prior to inflation would be inflated away |1Q . 



B. Quantum stress energy tensor 



We shall obtain some general results for the expectation value of the stress energy tensor using perturbation theory. 
The classical action of the radiation field 6 is 



S r = \ ( (j)A<j)dv (11) 
2 JM 



where A is a second order operator and dv is the volume measure. The unperturbed system is chosen to be in vacuum 
state where the green functions of the radiation with the prescribed boundary conditions are known explicitly. For 
scalar fields, the boundary conditions can be Dirichlct, where the field vanishes on the boun dary , or Robin, where the 
normal derivatives of the field are prescribed. Our notation for Green functions follows ref. |17| . 

The relation between the perturbed and unperturbed boundaries can be defined by a diffeomorphism / as shown 
in figure^ If the perturbation is small, then the vacuum state should be invariant. Expectation values calculated for 
the brane S and the metric g should be equivalent to the expectation values calculated with the unperturbed brane 
and the metric f*g. By this means we can replace a perturbation of the brane by a perturbation of the metric. 

The number of dimensions n is left arbitrary. Local coordinates can be adapted to the unperturbed brane so that 
n — 1 coordinates x^ describe the position on the brane and a coordinate z = —x n describes the distance to the 
perturbed brane. The perturbed brane is then described by the function z(x IJ '). 

We begin with the Wightman functions G ± of the scalar field <fi, defined by 

G+(x.x') = (m\(j>(x)(j>(x')\m) (12) 
G-{x.x') = (in\(f>(x')<f>(x)\m). (13) 

The Wightman functions satisfy the homogeneous wave equation 

AG* = (14) 

Their sum defines the Hadamard function, 

G (1) =G + + G~ (15) 
The change in Hadamard's function when the operator is perturbed by the moving brane is 

<5G (1) = -G fl (<5A)G (1) - G {1) {SA)G A (16) 



FIG. 1: The brane E can be obtained by mapping the unperturbed brane Eo along the intergral curves of a vector field £. The 
mapping is denoted by / and the pullback mapping by /*. 



where Gr and G a are the retarded and advanced propagators, 

G R (x,x') = G A {x',x) = -i(m\[<j>(x)<j>(x')}\in)6(t-t') (17) 

The perturbed form of the Hadamard function is real, symmetric and causal. 

The first order change in the Hadamard function can be cast into a form which is more convenient for calculations, 

SG {1) = -41m(G5 > (<5A)G^) (18) 

where the subscript denotes the use of the unperturbed green function, and 

G$(x,x') = G$(x',x) = iG+(x,x')6(t-t') (19) 

The subscripts will subsequently be dropped. 

The next step is to find an expression for SA. Consider the general second order operator 

A = -V 2 +eR + m 2 (20) 

The change in the operator is induced by the coordinate transformation £ a = — z{x^)n a as shown in figure ^ The 
diffcomorphism invariance of the classical action of the radiation 111|) implies that 



(</>(<5A)0 + (C a V a 0)A0 + <j)A(CV a 4>)) dv + / z<j)A(j)dv = (21) 
After integration by parts we deduce that 

SA = V a CA - ACV a (22) 
The required combination of green functions can now be found explicitly by a covariant volume integral, 

5G {1) =Alm( G > {x 1 x'){A£ > a '\/ a , -\/ a ,^ a ' A)G > {y 1 x')dv' (23) 

J M 

The imaginary part of the integrand consists of a total divergence which reduces by the divergence theorem to 

5G (1) =4Im f (S7 n/ G > (x,x')S/ n ,G > (y,x / )-G > {x,x')\7 2 n> G > {y,x / )) z{x')dv' (24) 



The reason for obtaining this expression for the Hadamard function is that it can be used to find the expectation 
value of the stress energy tensor. The expectation value of the stress energy tensor contains terms such as 

(in|(V a 0V b + V b 0V a 0)|in), (25) 
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which can be evaluated by applying an operator to the hadamard function and taking a coincidence limit |17|. 

[V V 6 ,G«]. (26) 
The brackets denote the coincidence limit x = x' . The full stress energy tensor is given by 

HTjin) = \[D a , b G^] (27) 



where the distribution valued operator D a b is given by 



1 S 2 „ 6A 
2\ 



Dab = o772^' ah = TTah' ( 28 ) 



In the minimal case, for example, 

1 



A,&(*> x 1 ) = (^V a , V b - ^9ab9 cd V c> V d j 5(y, x)S(y, x') (29) 

For Dirichlet boundary conditions, the change in the normal component of the stress energy tensor at the brane 
position is given by (|24p. 

5<in|T^ n |m) = Im / (V n V n >G>(x, x')) 2 z{x') dv' (30) 



which is equivalent to an old result due to Ford and Vilenkin [(j. For a moving mirror, this stress energy component 
generates the pressure force on the mirror. Since the unperturbed value of the stress energy does not affect the motion 
we shall drop the S and only keep the perturbed value from now on. 



C. Schwinger-Keldeysh version 

Before proceeding, it is interesting to see how the same results can be obtained more directly from an effective 
action using Schwinger-Keldeysh methods [l8L [Tgt |20| . In this approach, the operator A is a 2 x 2 matrix which 
depends on two copies of the background fields, 



1SK 



AM o 

A[g 2 ] 



(31) 



The radiation reaction forces can be obtained by variation of the one loop correction to the Schwinger-Keldeysh 
effective action, 

W SK = - l -logdet{iG SK ) (32) 

where Gsk is the green function in the Schwinger-Keldeysh approach. Let S[g±] be the classical gravitational action, 
then S±S + 8\Wsk — Q> where 8\ denotes variation of g±, and g = g\ = g 2 after the variation. 
The Schwinger-Keldeysh green function is a 2 x 2 matrix 

GsK= ( G zG+ -Gf ) (33) 

where 

G T (x,x') = iG + 9{t-t')+iG-0(t' -t) (34) 

G T {x.x') = iG-9{t-t') + iG + 6(t' -t). (35) 

(Note that Gt is identical to the Feynman Green function if the in and out vacua are the same state, but not 
necessarily identical otherwise) . The Schwinger-Keldeysh Green function changes under perturbations of the operator 
by 



&Gsk — —Gsk{SA)skGsk 



(36) 
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Multiplying out the matrices gives, for example, 

SG T = -G T (6A)G T + G~ {8A)G+ (37) 
Under the same restrictions as before, this is equivalent to 

SG T = -2 Im(G > (5A)G< ) (38) 
The variation of the Schwinger-Keldeysh effective action to first order in perturbation theory is therefore 

SiWsk = lmtr(((5iA)G > ((5A)G < ). (39) 

By (|27|l and l|28|) . this is equivalent to the previous result for the stress energy tensor of the radiation. 

Alternatively, we we can regard the action as a function of the brane position and obtain a reduced action formu- 
lation, often called the moduli space approximation. Let 

W[z,g] = -Imtr (5AG > 6AG<) . (40) 

Variation with respect to z can be related to a metric variation 5{W by diffeomorphism invariance, as in figure ^ 
The effective equations of motion are then 

oz oz 

The Schwinger-Keldeysh formalism is usually applied to non-equilibrium thermal field theory, using finite tem- 
perature version of the propagator (|33|) . Our results are therefore equally applicable to finite temperatures (for a 
time-independent background) if we replace the function G > {x 1 x') by the finite temperature version. 

III. BRANES MOVING IN FLAT SPACE 

We shall examine the radiation back reaction force on the brane to leading order in the displacement from a 
fiat hyperplane in flat space. Consider the massless field with Dirichlet boundary conditions to begin with. The 
unperturbed green function l|19|l can be expressed in terms of basis functions which vanish at z — 0, 

G>(x,x') = jT J ^^^Asinqzsinqz'e^-^-^-^eit-t') (42) 

where u> — (k 2 + q 2 ) 1 / 2 . Inserting this into equation (|3(J|) for the stress energy tensor gives 

m T„ n m) =Imreg / — — / r I[z] (43) 

J_ 00 2tt 2tt J (2tt)" 1 ujuj' 

where 'reg' indicates some form of regularisation has been performed to make the integral finite and 

I[z}= f e-^+^-t'h^dt' = -i((w + w')-id t y 1 z (44) 

J —OO 

Note the importance of taking the imaginary part of the expression after regularisation. It is useful to introduce 

f dg dq> d n - 2 k q 2 q 12 
n[X) ~ r6g J 2tt 2ir (2tt)"- 2 u>lj' (w + w') - ix 1 ' 

and then 

(m\T; m \m)=lmF n (d t )z (46) 

The integral diverges for all values of the dimension n which rules out dimensional regularisation. However, an 
analytic regularisation scheme can be used where we define 

F (x) 1 r dXy f ^ dn " k (47 ) 

* ns[X) ~ T( S + l)J AA J 2tt 2ir (27t)»- 2 uu> 1 > 
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and take the value at s = 0, removing pole terms if necessary. After integrating over q, we have 

F ^ = r^TT)l ^7(^ (Afer ^ l(Afc)2e (48) 

where K\ is a Bessel function of the second kind. The remaining integrals give 

F (x) - 1 *"" S+1 F(n/2)3 (49) 
nsK ' r(s + l) 2W 2 (n + l)T{n) 2 smTT(n-s) K ' 

After removing the pole at s = 0, the regularised expression is 

(-ix) n+1 ( , (x\\ T(n/2) 3 . . 

where /i is a renormalisation constant. 
Some examples are 

ImF 2 (x) = (51) 
*n*i(*) = ^ (52) 



16384tt 2 "* V M 



^(z) = -— —-ln(^) (53) 



The results for two and four dimensions agree with those found by Ford and Vilenkin for the moving mirror problemlyl . 

The dependence on /i in odd dimensions can be traced to need for new counterterms in the classical action |2l| 
which arise from the n'th heat kernel coefficient of the operator A j^. In five dimensions, the heat kernel coefficient 
contains boundary terms of the form 

S 5s = - [ V^K^V^K^dv, (54) 
s Js 

and other permutations of the indices, where K is the extrinsic curvature of the brane. These counterterms produce 
a dt 6 z in the equation of motion (since K = ~dt 2 z to leading order) which cancels the divergence in the radiation 
reaction and leaves the logarithmic dependence on /i. How the logarithmic derivatives are interpreted, and how to 
solve the equations of motion for the brane will be covered in section 4. 

The radiation back reaction from a massive scalar field with mass m can be found by redefining uj = (fc 2 + q 2 +m 2 ) 1 ^ 2 
in eq. (|47|l . In even dimensions, the integral results in 

ImF n (x)= - \ , m"~Vrf— Vf-Wi f— ,-;2:-- ^V) (55) 
V ' 8(47r)"/ 2 V 2 / V 2 / V 2 2 4m J 



where 2-F1 is a hypergeometric function. In odd dimensions, we discard the pole term and retain only the finite part 
of the same expression. 

For small values of x, F n is of order x 2 . If the position of the brane is oscillating with a frequency u> < m, the 
radiation reaction force is proportional to ui 2 , which is similar to a known result for a domain wall in four dimensions 
|23|. For large values of x, the hypergeometric functions have branch cuts in the region |x| > 2m and logarithmic 
terms appear in the large x limits for both even and odd dimensions. The leading terms agree with the massless 
results given above. 



IV. BRANES MOVING IN ANTI-DE SITTER SPACE 



Branes in anti-de Sitter space are interesting from a cosmological point of view. The intrinsic geometry of a 
homogeneous brane moving in anti-de Sitter space is similar to a cosmological model 0, |2^. The effects of the 
higher dimensional cosmological constant and the vacuum energy on the brane can be fine-tuned to give a relatively 
small effective cosmological constant on the brane [2^, 12^ ■ 

As before, we will neglect the back reaction of the motion of the brane on the bulk metric when calculating the 
quantum stress energy tensor. We shall consider small perturbations of a single flat brane whose penrose diagram is 
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FIG. 2: On the left is a Penrose diagram of universal anti-de Sitter space showing the location of a timelike slice of the flat 
brane. The coordinate system extends on the right of the brane (z > a) to the horizons Tt + and HT where the timelike killing 
vector becomes null. On the right is a copy of the Penrose diagram to show the location of a de Sitter space brane. 



shown in figure [3 We take the vacuum state related to the timelike translation symmetry along the brane, which we 
call the brane vacuum. 

The brane vacuum has the disadvantage that there are horizons where the timelike killing field vanishes. The 
horizons can be avoided, for example by considering a two brane system and taking the limit where the separation 
between the branes becomes large. If, instead, we use a vacuum state which covers the whole of anti-de Sitter space, 
then we have to consider the effects of Bogolubov coefficients ^t|- I n the case of the de Sitter brane, the Euclidean 
vacuum state can be used |2l| . The Euclidean vacuum becomes a thermal state with respect to the brane, creating 
additional thermal effects in the radiation reaction force. We shall not consider this further here. 

One of the many ways to express the Anti-de Sitter metric in n dimensions is the conformally flat form, 



a z 



~ 2 (dz 2 + rj f , u dx fi dx v ) (56) 
with Minkowski metric r]^. The Anti-de Sitter radius a is related to the cosmological constant in n dimensions by 

(n-l)(rc-2) 

A " = 2^ • (5?) 

The unperturbed brane will be placed at z — a and the manifold extends in the z > a direction only. The covering 
space can be filled in by reflection symmetry about the brane. 

The Green functions can be obtained from the normalised modes Uk q of the scalar wave equation. These are 

^=(^g) ( "" )/2 ^(^K (k " } (58) 

with u> = (k 2 + q 2 ) 1 / 2 and H v a combination of Hankel functions of order 

1 



p 2 



A +a 2 {(t-i(n))R + m 2 ). (59) 



Conformal curvature coupling corresponds to the case £ = £(n) and v = 1/2. 
For modes which vanish on the brane, 

H v (q,z) = (H^( q z)Hi 2 \qa) ffW (qz)H^ (qa)) (60) 

Note that the normal derivatives of the mode functions when evaluated on the brane z = a are then given by 

V„«,= ' (^) " 2 -^3- e«— > (61) 
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In the conformal case v = 1/2, this reduces to 

v n^kq 



(2w) 1 /2 



2 1 /2< ? gi(k-x-wt) 



(62) 



which is also the value obtained in flat space. 

The normal derivatives of the Wightman function are needed to obtain the stress energy tensor. The Wightman 
function can be obtained from a mode sum, 



dq f d n -^k 



2tt 



G + (x,x') -- 

The normal derivatives at the brane can be expressed as 

dq f d n - 2 k 2q 2 
)«- 2 2uj 



(2tt)"- 2 WzK 9 0*O 



(63) 



H[%(qa) 



H^\qa) 



i(k-(x-x')-w(t-t')) 



(64) 



As in the flat space example, the expectation value of the normal components of the stress energy tensor are given in 
terms of an integral by 



(in|r„ r Jin) = lmF n (d t )z 



where the anti-De Sitter space version of F is 

dqdq' d n ~ 2 k q 2 q' 2 



F n (x) = rcg 



2tt 2tt (2tt)"- 2 ulj' 



H^ 2 (qa)H^ 2 (q'a) 



H^\qa)H^>{q'G 



(uj + lu') — ix 



(65) 



(66) 



The first striking feature of eq. I|66|l is that, in the conformal case v = 1/2, the result for radiative back reaction is 
identical to the flat space result. In the non-conformal case, we can obtain the large ax limit from the large argument 
expansion of the Hankel functions, for example 



ImF 5 (z) 



In I* 



163847T 2 



(67) 



in five dimensions. The small ax limit can correspondingly be obtained from the small argument expansion of the 
Hankel functions, which leads to 



ImF 5 (x) - Cx b (ax) 



(68) 



when ax <C 1, where C is a constant. 



V. EQUATIONS OF MOTION WITH LOGARITHMIC TERMS 

We can now construct the equation of motion for the brane. For a small perturbation z of a flat brane, the scale 
factor awl and d t z « —aH. For conformal scalar field radiation, eq. l|6T|l suggests that eq. © would become 

-dF + 2H - 2H ^- AK ^{-,dt)^- (69) 

where Hq = A/3 and A is a numerical coefficient. Strictly speaking, this equation is incomplete because there may 
be additional radiation damping terms of order H 2 . 

The radiation damping term calculated here is small in recent cosmological eras. Returning to eq @, we see that 
the combination n 2 a ~ Ga 2 in order of magnitude. The experimental lower bound on (Ga 2 )^ 1 / 4 is around ITeV, 
but there is no reason for the the value to be much smaller than the usual Planck scale. When H ~ H 2 , the damping 
term is only important when H 4 ~ (Ga 2 )^ 1 . 
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There is a small residual dark radiation effect given by eq. I|10|) . If we take into account a period of inflation, ending 
at time i/, then the ratio of dark radiation to ordinary radiation is 

n a 2 f°° H , fld\d 5 H J 

6AkI<j / — In — - —nrdt 70) 
5 J tl p \fidtj dt* V 1 

which is of order Ga 2 H(tj) 4 . In principle, H(tj) should be close to the value of H determined by the tensor 
perturbation amplitude Ax in the cosmic microwave background, A T = 8nGH 2 |28|| . To be within nucleosynthesis 
constraints, the amount of dark radiation must be small fl2. Il6|. and therefore there is a lower bound on the string 
vacuum energy scale of approximately 

(Ga 2 r 1/4 > A~ 1/2 m p (71) 

where m p is the Planck mass. 

We shall analyse a slightly more general type of equation of motion, 

dp , , / 1 d \ d n p , , 

-i--F(p) = -T]ln[-—)—^ 72 
dt w ' \pdtj dt n ' 

A reasonable definition of the logarithmic term should be linear and causal, to ensure that the back reaction depends 
on the history of the source. We can begin be defining the action of the logarithmic derivative on exponentials, 



1 d 



l "{ji t )<°'= l »{-»)'°' (73) 

where we take a > 0. This leads to the following definition, 

hi (^)^-£ i "< e, "< , - t '»l < «' < 74 » 

where 7 is Euler's constant. The definition is linear, causal and correctly reproduces the action on exponentials. The 
integral can be evaluated for a broad class of functions, including discontinuous functions, for example 



1 d 



a 



In y-j t J 6(t)e at = I Ei (erf) +ln l^-jj 6{t)e at (75) 

where Ei is the exponential integral. If / is integrable, then the logarithmic derivative of / is 0(t~ r ) as t — > 00. 
For a linear equation of motion l|72|) . with F — — Xp, we have 

The exponential solutions exp(pzt) satisfy 

z + a = -I3z n ln(z) (77) 

where a — p~ 1 \ and j3 = ■qp: n ~ 1 . The roots for a sample case are plotted in figure|3 When there is no reaction term 
there is only a decaying mode corresponding to z = —a. The reaction term shifts this solution and introduces an 
oscillation. There are two other decaying solutions and two growing modes. These modes have \z\ > 1 and damping or 
growth rates larger than the renormalisation scale p. We tentatively identify these as unphysical runaway solutions. 

Runaway solutions can be excluded by replacing the equation of motion by an integrodifferential equation. First, 
we introduce a Green function which decays exponentially as the time t — > ±00 and satisfies the following equation, 

(1 rl \ rl n ~ 1 r* 



The equation of motion (|72|l can then be rewritten in an alternative form 

G(t-t')F{t'). (79) 



dp 
~dl 
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FIG. 3: The roots of 1 + z = — 0.1z J ln(z) in the complex plane. 



The runaway solutions are excluded at the expense of a preacceleration term, as we can see from the expression for 
the Green function 

G(t-t')= ^ Cie^-^eit-t')- Cie^-^Otf-t), (80) 

Rezi<0 Rezj>0 

where Zi are the roots of Tqz n ~ l ln(z//i) = — 1 and Cj are constants fixed by continuity relations at t = t'. The 
acausality is represented by the final term, but it is a relatively minor effect if Rez^ 3> /i, because then the acausality 
only occurs on timescales small compared to the renormalisation timescale. 

An interesting variant of the equation of motion l|72|l occurs for the case n = 1 and F = 0. There are two real 
solutions, p — and p = exp(exp(— Because of the small size of the exponent it would not be appropriate to 
regard the exponential solution as a runaway solution. In a cosmological context, a flat brane would spontaneously 
begin to accelerate with a tiny, positive cosmological constant. We have not been able to find a model with this 
particular radiation reaction so far. 



VI. CONCLUSION 



This paper has been devoted to calculating the radiation back reaction forces on a moving brane. We have found 
that massless, conformally invariant scalar radiation into the bulk results in logarithmic terms in the equations of 
motion. The radiation reaction leads to both physical and runaway solutions. The runaway solutions can be excluded 
by introducing an integrodifferential formulation of the equation of motion. This form of the equation could be used 
if we were to attempt to solve the equations of motion numerically. The question of whether the brane's initial 
trajectory might be destabilised by radiation reaction forces remains an open one. 

The radiation from the moving brane also leaves a residual effect in the form of dark radiation. The amount 
of this radiation depends on the brane trajectory |J7DJ|. An order of magnitude estimate combined with a crude 
nucleosynthesis constaint gives a lower bound on the string vacuum energy scale of V 'At m p , where At is the tensor 
perturbation amplitude in the cosmic microwave background. 

The radiation reaction problem has been set up in a way which can be generalised to non-scalar fields. One issue 
which would be of interest is the possibility of supersymetric cancellations amongst some of the terms. The regulari- 
sation methods used in this paper have enabled us to neglect the non-logarithmic divergences, but it would be better 
to see direct cancellation of these divergent terms. For broken supersymmetry, we expect that the renormalisation 
scale, which appears in our results, would be replaced by a supersymmetry breaking scale. 

The results for the radiation reaction of branes moving in anti-dc Sitter space generalises in a very staightforward 
way to any background where the modes of the radiation fields are known explicitly. An example of this is the low 
energy limit of the heterotic string 0, |2^, [3(j , where the modes of the graviton multiplet are known and can be found 
in [p. 

Another aspect of the radiation reaction which is worth further study is the effect of different higher dimensional 
vacuum states. These states can have thermal properties and the reaction forces would have similarities to those 
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acting on moving mirrors at finite temperatures. 
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